Nematic superconductivity in magic-angle twisted bilayer graphene from atomistic modeling

Twisted bilayer graphene (TBG) develops large moir\'e patterns at small twist angles with flat energy bands hosting domes of superconductivity. The system size and intricate band structure have however hampered investigations into the superconducting state. Here, using full-scale atomistic modelling local electronic interactions, we find above experimentally relevant temperatures a highly inhomogeneous state nematic ordering on both atomic length scales. has locally anisotropic real-valued d-wave pairing, vector winding throughout pattern, is three-fold degenerate. Although symmetric, full gap, which tie to {\pi}-phase interlayer coupling. nematicity further directly detectable in density states. Our results show that modeling essential also very similar interactions produce different states TBG high-temperature cuprate superconductors.